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Abstract — The nonlinear surface impedance of a thin su- 
perconducting strip carrying a microwave current has been 
calculated numerically from first principles based upon flux 
penetration due to a current-induced critical state. The cal- 
culations approach analytic results in the appropriate limit. 
The technique is useful for design of passive superconduct- 
ing microwave circuits. 
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I. Introduction 

SINCE its discovery, the new breed of oxide supercon- 
ductors with transition temperatures above the boiling 
point of nitrogen have shown great potential as a replace- 
ment for metals for passive microwave structures owing to 
their very low surface resistance^, 0, |J. However, a key 
feature of the HTSC is the inherently nonlinear response at 
microwave frequencies, which is present over a wide range 
of input powers. This nonlinear responses also poses a ma- 
jor challenge for numerical device design using these ma- 
terials. Most existing computer aided design software do 
not take into account the inherent nonlinearity of super- 
conductors. The few analytical models that do exist B are 
valid only for highly idealized geometries and need to be 
incorporated into practical device design. 

There are several physical models of nonlinear surface re- 
sistance, each valid in a different regime of temperature and 
microwave power The most pursued ones are pinning 
induced critical state Q, S, 0, transport current induced 
self- heating || , and depairing B . The current-driven crit- 
ical state model has proven to be particularly useful not 
only because it fits experimental measurement of Rs very 
accurately over a wide range of temperature and power lev- 
els ||, but also because its usefulness extends to modeling 
hysteretic losses in superconducting transmission lines. 

The (field-induced) critical state model was first pro- 
posed by Bean |]lo|| to explain the unusual magnetization 
behavior of hard superconductors and was first solved for 
the geometry of a semi-infinite slab where the field was ap- 
plied parallel to its face. The geometry where the field is 
applied perpendicular to a thin strip was solved relatively 
recently by several groups and the field and current distri- 
butionjTlJ, E^] and the surface impedance from hysteretic 
lossQ was obtained analytically. 

Though very elegant, the analytic results were restricted 
to certain "ideal" geometries only where symmetry could 
be used to simplify the problems. However, these idealiza- 
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Fig. 1. The geometry of the problem. The strip is centered around 
x = y = 0. 



tions are not always satisfied in real world applications of 
superconducting strips in microwave circuits. In the wake 
of the fact that the model does describe experimental data 
accurately 0, it is imperative to be able to implement a 
numerical paradigm without the constraints of the ideal- 
ized geometry required by the analytic approach. 

In this paper, we demonstrate a numerical approach to 
incorporate the nonlinear response of superconducting ma- 
terials using assumptions equivalent to those of the current 
induced critical state. 

II. The Critical State Model 

Let us consider a striplinc with the ground plane far 
away from the center conductor so that it carries a negli- 
gible current and hence does not effect the magnetic field 
distribution inside the center conductor. It has been shown 
in literature that such an assumption is reasonable [|l3| (We 
will discuss the effects of ground planes somewhere else.) 
The geometry of the center conductor and our choice of 
axes is shown in of Fig. 1. 

Starting from the most generic physical assumptions that 
retain the essential ideas of critical states: 

1- |j( r )| < jc where j is the current density. 

2. The magnetic induction B(x, y) = for \y\ < b < a 
The first assumption highlights the fact that flux lines are 
pinned until the Lorenz force on them exceeds a certain 
critical value. The second assumption follows from the fact 
that the vortices start penetrating from the edges. Hence 



for any given current there would be a flux free region in 
the center of the strip. Note that the depinning critical 
current is at least an order of magnitude smaller than the 
depairing critical current where superconductivity would 
be quenched by the self field. 

In the following we show that for a strip with finite thick- 
ness, this model can be solved numerically and that it leads 
to unique field and current distribution within the strip. 
We further show that (1) and (2) leads to hysteresis irre- 
spective of the geometry and study the hysteresis laws for 
the specific strip geometry under consideration. 

III. Computational Details 

A semi-numerical approach was used to determine a cur- 
rent distribution which satisfies the two assumptions pre- 
sented above. Assuming an uniform current distribution 
along the x-direction one can show that 




Fig. 2. Penetration depth b as a function of current I //max- 
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The above analytical expression for B(x — 0, y) in Eq(Q) 
was used to save computational time. 

In the model which requires that j(y : b < \y\ < a) = J c 
and B(x = 0, \y\ < b < a) = 0, one needs to solve for the 
current distribution in the field free region \y\ < b < a 



j(y')tan~ 



2(y -y') 



-]dy' = for \y\ < b < a (2) 



Numerical solutions of integral equations such as eq(|2|) 
are rather complicated. We found a simpler iterative ap- 
proach to find the current density in the field free region of 
the superconducting strip. Here we give a brief description 
of the algorithm that was used. 

Suppose that the thickness of the strip is infinite (5 — > 
oo.). Then for any given initial current distribution j(y), 
a field free region B(\ y\ < b) = can be obtained by 
replacing the current density in that region by 



m<b)^j{\y\) + -^B x 



(3) 



Also note that the magnetic field strength on the line x = 
for a slab (5 — > oo.) is mostly effected by the current on 
and around the line x — 0. Therefore, for the case of a 
finite 6, the transformation in eq([|) will bring the field 
strength in the region \y\ < b a little closer to zero. Hence 
an iterative implementation of the above transformation of 
J(y) in eq(^|) will finally converge bringing the magnetic 
field strength in the region \y\ < b to zero, irrespective of 
the size of 8. 

In a nutshell, we start with a flat current distribution, 
j(y) = j c in an one dimensional grid of size 4000 and we 
iteratively make B(x = 0, \y\ < b) vanish by replacing the 



current density by j™ 
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IV. Results 




tower (darker) -Finite thickness 
' (Numerical calc.) 
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Fig. 3. Current distribution j/jc across the strip for b/a = 
0.1, 0.2, . . . , 0.9. The lower (solid) lines are numerical calcula- 
tions and the upper (dashed) lines are analytic calculations. 



A. Penetration Law 

We calculate the flux penetration depth d/a for a strip of 
aspect ratio of 1% as a function of the applied current. In 
Fig. 2 we display (1 — d/a) = 6/a as a function of the applied 
current y/a and compare that with the analytical results 
for a strip of zero thickness (T^]. Fig. 3 and Fig. 4 shows 
the corresponding current and field distributions. One can 
see from Fig. 2 and Fig. 3 that the penetration law get's 
more linear and the current density in the field free region 
decreases as we go from zero thickness to a finite thickness. 
These results are expected since, in the limit S — > oo we 
expect a linear penetration law and J{x — 0, \y\ < b < a) = 
0. These results can then determine the current critical 
states and the corresponding magnetic field distribution 
for AC currents with a peak current of Iq. 

V. Hysteresis 

The hysteresis law can be derived by studying the current 
density distribution for an oscillating current over a full 
cycle. If the current density distribution for a total current 
/ is j(y,I,j c ~^ °°) then in absence of pinning the current 
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Fig. 4. The magnetic field strength B(x = 0,y/a) across the strip 
for /// m ax = 0.41, 0.58, 0.70, 0.79, 0.86, 0.91, 0.95, 0.98, 0.99. 



density will be Nj(y, I,j c 
In other words 



ool for a total current of NI. 



j(y,NI,j c -> oo) = Nj(y,I,j c -> oo) 



(4) 



However, since j(y) < j c everywhere, the current density 
distribution depends on the total current. That is 



j(y,NI,j c )^Nj(y,I,j c ). 



(5) 



When the current is reversed, after reaching the peak I , 
current density starts to decrease from the edges with an 
effective critical current density in the opposite direction 
of 2 J c . This is because at the edges the current density is 
j c in the positive x-direction. 



3d = 2j 



(6) 



Thus the hysteresis law on the downward cycle is given 

byEi 

3l(v>I)3c) = jmax(y,Io,jc) ~ 3^{yJo ~I^3c) (7) 

Fig. 5 displays the magnetic field distribution when the in- 
put current is reduced from the maximum current 2J c Sa to 
—2J c 5a. As the decreases through one does not see any 
field free region inside the strip, that is there are locked 
in fields from field distribution when the total current was 
2J c da: something that is unique to hysteresis. 

VI. Surface Impedance 

Assuming that the strip is excited with a current of am- 
plitude io and frequency ljq, the electric field E(y,t) is be 
given by 



(8) 



E z (y,t) = - J dy' B x (y',t) 
and the energy lost in a half cycle is 



U H .i 



2 / dy 
'0 



T/4 



T/4 



dt J(y,t)E(y,t) 



(9) 
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Fig. 5. 

from /max = 2J c Sa to —/max- This graph dispalys the ori- 
gin of hysteresis. We plot B across the strip for ///max = 
0.18, -0.17, -0.40, -0.58, -0.72, -0.82, -0.90, -0.96, -0.99. 



Substituting eq(ph and integrating by parts we get 



U H ., 



P T ' i d 
2 I dy dt J, 



-T/4 



dt 



dy' B x {y' ,t\\0) 



b(I ) 



pa pa 

4/i J c / dy dy'B x (y',t) (11) 

Jb(I ) Jb(I ) 



This definition of Ur.c. wa s also used by Brandt in ref . p2| . 
This defines the surface resistance as 

R s = 2^ P dy f dy'B x (y',t) (12) 

1 Q 1 Jb(Io) Jb(I ) 

Fig. 6 displays the surface resistance R s as a function of the 
peak current, Jo- 

The surface reactance X s per unit length can be calcu- 
lated using the average total electrostatic energy per unit 
length. 



A, = 



(13) 



We carry out the space integral until B 2 (x, y) falls down 
by 99% of it's peak value. Fig. 7, displays the surface reac- 
tance X s as a function of the peak current, Iq. Our results 
are close to the analytical results in ref.|Q. 

VII. Conclusion 

A first principle numerical calculation on a strip of fi- 
nite thickness using assumptions similar to the critical state 
model yields results that approach the analytic results in 
the appropriate limits (6 — > 0) but without the unrealistic 
idealizations of the later. Further a numerical approach 
gives one the flexibility of accommodating higher complex- 
ities like effects due to the ground planes and more ac- 
curate calculations of surface impedance and harmonics. 
This topic will be addressed in future work. This approach 
can be implemented into a computer aided design frame- 
work for the design of passive superconducting microwave 
components. 



4 



[8] M. Golosovsky and M. Tsindlekht and H. Chayet and D. Davi- 
dov "Vortex depinning frequency in YBa2Cu307- x supercon- 
ducting thin films — anisotropy and temperature dependence," 
Phys. Rev. B, vol. 50, no. 1, pp. 470, Jul 1994. 

[9] Balam A. Willemsen, John S. Derov, Jose H. Silva, and S. Srid- 
har, "Nonlinear response of suspended high temperature super- 
conducting thin film microwave resonators," IEEE Trans. Appl. 
Supercond., vol. 5, no. 2, pp. 1753-1755, June 1995. 

[10] C. P. Bean, "Magnetization of hard superconductors," Phys. 
Rev. Lett, vol. 8, pp. 250-253, 1962. 

[11] E. Zeldov, John R. Clem, M. M. McElfresh, and M. Darwin, 
"Magnetization and transport currents in thin superconducing 
films," Phys. Rev. B, vol. 49, no. 14, pp. 9802-9822, April 1994. 

[12] Ernst Helmut Brandt and Mikhail Indenbom, "Type-II- 
superconductor strip with current in a perpendicular magnetic 
field," Phys. Rev. B, vol. 48, pp. 12893-12906, 1993. 

[13] A. M. Campbell, " " IEEE Trans. Appl. Supercond., vol. 5, pp. 
687, 1995. 



Fig. 6. Surface resistance R s as a function of the current ///max- 
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Fig. 7. Surface reactance X s as a function of the current I/Imax- 
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